clear all
%load ParEst_eh0_time0_Mex_nokids ParEst_eh0_time0_Mex_nokids
load ParEst_eh0_time1_Mex_nokids ParEst_eh0_time1_Mex_nokids
%Chute = [1.5,1.5,ParEst_eh0_time0_Mex_nokids(17:20)];% % gamma_h, gamma_l, alphaf_2, betaf_2, alphai_2, betai_2
Chute = ParEst_eh0_time1_Mex_nokids;
HandleLo = @Principal_A2;
LoBo = [-1, -1, 0.1, 0.1, 0.1, 0.1];
HiBo = [2.5, 2.5, 10, 10, 10, 10];
%options = optimoptions(@simulannealbnd, ...
 %                    'PlotFcn',{@saplotbestf,@saplottemperature,@saplotf,@saplotbestx,@saplotx},...
  %                    'FunctionTolerance',10,'ReannealInterval',2, ...
   %                  'InitialTemperature', 2*10^4*[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 10 10 10 10 10 60], ...
    %                 'TemperatureFcn', @temperaturefast, 'MaxStallIterations', 4000);
                 options = saoptimset( ...
                     'TolFun',3000,'ReannealInterval',4, ...
                     'InitialTemperature',10^2*[3.5 3.5 10 10 10 10], ...
                     'TemperatureFcn', @temperaturefast, 'TimeLimit', 2*86400); % temp default 10*10^2* 
[x,fval,exitflag,output]  = simulannealbnd(HandleLo, Chute, LoBo, HiBo, options)
ParEst_eh0_time1_Mex_nokids = x;
save ParEst_eh0_time1_Mex_nokids ParEst_eh0_time1_Mex_nokids